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INTRODUCTION 


When  analyzing  engagements  between  submarines  and  surface  ships,  it  is  often 
necessary  to  estimate  the  effectiveness  of  the  submarine's  torpedo  attack.  Although 
sophisticated  homing  torpedoes  are  being  used  by  many  of  the  world's  navies,  there 
continue  to  be  large  inventories  of  straight-running  torpedoes.  This  research  con- 
tribution presents  a methodology  and  a set  of  computer  programs  for  analyzing  an 
attack  on  a nogpvading  surface  ship  by  submarine(s)  firing  straight-running  torpedoes. 

It  explicitly  assumes  a triangular  horizontal  cross  section  for  target  ship  geometry. 

We  believe  this  approximation  is  more  applicable  to  combatants  than  previously  assumed 
geometries  and  provides  reasonably  accurate  results  for  bow  and  stern  shots  as  well 
as  beam  shots.  The  programs  provided  make  it  possible  to  analyze  the  effects  of  several 
simultaneous,  independent  submarine  attacks.  The  user  must  provide  some  of  the 
torpedo  and  target  ship  characteristics  and  submarine  firing  doctrine.  These  are: 


Torpedo 


- speed  (knots) 

- maximum  run  distance  (yards) 

- reliability 


Target 


- length  at  waterline  (feet) 

- maximum  beam  (feet) 

- speed  (knots) 


• Firing  doctrine 


- number  of  torpedoes  in  the  salvo 

- salvo  spread  (expressed  as  a multiple  of  ship  length). 


The  submarine's  position  relative  to  the  target  is  a control  variable.  The  position  is 
expressed  as  a range  and  bearing  measured  clockwise  relative  to  the  target's  heading. 


If  a vulnerability  curve  is  available,  it  is  possible  to  compute  the  probability  that 
the  target  sustains  a particular  level  of  damage.  The  David  Taylor  Ship  Research  and 
Development  Center,  Carderock,  is  one  source  that  provides  vulnerability  curves  for 
a ship  being  put  out  of  action,  sunk,  or  immobilized. 


MODELING  A TORPEDO  ATTACK 


The  Torpedo  Problem 

The  classical  torpedo  problem  is  in  three  parts: 


1- 


(1)  Given  that  a salvo  of  n torpedoes  is  fired  at  a ship,  calculate  the  probability 
that  0, 1, 2, . . . , n torpedo  tracks  intersect  the  target  if  all  torpedoes  function 
properly. 

(2)  Given  the  results  in  (1),  calculate  the  probability  that  0, 1, 2. . . , n torpedoes 
reach  the  target  and  explode. 

4 , 

(3)  Given  that  0, 1, 2, . . . , n torpedoes  hit  the  target  and  explode,  calculate  the 
probability  that  the  target  is  out  of  action,  severely  mission  degraded,  or  sunk. 

To  solve  the  first  part,  consider  the  typical  fire  control  triangle  shown  in  figure  1. 

It  is  easily  seen  from  this  figure  that  torpedo  lead  angle,  X , is 

X = arc  sin  sin  pl  . (1) 


FIG.  1:  TYPICAL  FIRE  CONTROL  TRIANGLE 


Although  x does  not  depend  directly  on  range,  errors  in  measuring  range  ( to 
estimate  u and  0 ) result  in  lead  angle  errors.  Lead  angle  error  is  derived  as: 


E = standard  deviation  of  the  lead  angle  error, 

K 

Eu  - standard  deviation  of  the  errors  in  estimating  u 
Ep  = standard  deviation  of  the  errors  in  estimating  j3 
E = standard  deviation  of  incidental  errors. 


Values  for  these  standard  deviations  have  been  estimated  in  reference  1 from  short 
range  (less  than  4,000  yards)  attacks  on  Japanese  ships  during  World  War  II  as  follows: 


10  degrees  (0.  1745  radians),  and 


E.  = 0.05,  0.03,  and  0.015  radians  for  1,000,  2,000,  and  3,000-y"ird 
ranges  to  the  target  at  the  time  of  firing. 

Other  values  of  E , where  needed,  are  taken  from  the  continuous  function: 


These  parameters  are  used  in  the  computer  program  presented  later,  Data  from 
more  recent  experience  with  straight-running  torpedoes  can  be  easily  inserted  into  the 
program,  if  it  is  available. 


The  Model 

E^  is  used  to  estimate  torpedo  hit  probabilities,  assuming  that  the  torpedo  track  has 

a Gaussian  (normal)  distribution  about  its  intended  course  with  mean  zero  and  standard 
deviation  E ^ . The  probability  of  a hit  on  a ship  by  a single  torpedo  aimed  at  the 

center  of  the  target  is  the  integral  from  -Leff/2E^RT  to  +Lgff/2E^RT  of  the  standard 

normal  distribution.  is  the  effective  length  of  the  target  measured  perpendicular 

to  the  torpedo  track,  and  RT  is  the  length  of  the  torpedo  track. 


For  a salvo  of  torpedoes,  the  limits  of  integration  for  0, 1, 2,  or  more  hits  are  a 
function  of  L^^/E^RT  , the  spread  of  the  salvo,  and  the  number  of  torpedoes  in  the 

salvo.  When  a salvo  of  more  than  one  torpedo  is  fired  at  a target,  it  will  be  assumed 
that  the  center  of  the  salvo  is  fired  at  the  center  of  the  ship,  the  torpedoes  are  equally 
spaced,  and  the  individual  torpedoes  within  the  salvo  run  in  parallel  with  no  error  be- 
tween torpedoes.  The  only  error,  E , is  the  placement  of  the  center  of  the  salvo  on 

the  target. 


A formula  for  E was  presented  in  equation  2.  From  figure  1,  it  is  apparent  that 


RT  = vt  . We  next  need  to  determine  L 


eff  * 


To  obtain  an  exact  formula  for  requires  a description  of  the  horizontal  cross 


section  of  the  target  ship  at  torpedo  running  depth.  This  is  impractical.  Prior  models 
generally  use  approximations  based  only  on  length,  with  the  beam  assumed  to  be  zero. 
This  produces  severe  errors  for  bow  and  stem  shots.  Often  implicitly  rather  than 
explicitly,  some  geometric  shape  is  assumed  for  the  waterline  cross  section.  In  many 
cases,  the  shape  is  a rectangle  of  length  L and  width  B,  clearly  applicable  to  barges 
and  to  a lesser  extent  merchant  vessels  and  auxiliaries,  but  not  to  combatants.  Modem 
combatant  hulls  tend  to  have  fine  lines  at  the  bow  and  nearly  uniform  width  aft  of  the 
maximum  beam. 


We  considered  several  geometric  shapes.  For  example,  the  combination  of  a rectangle 
plus  an  isosceles  triangle  provides  an  excellent  approximation,  but  a combination  of 
geometric  shapes  complicates  the  mathematical  formulation  considerably.  In  general, 
compound  shapes  result  in  complex  mathematics.  A triangle  alone,  however,  provides 
results  that  are  nearly  as  good  with  much  simpler  mathematics.  Therefore,  we  used  an 
approximation  based  on  a triangular  shape,  with  the  small  angle  pointed  forward  and  the 
base  as  wide  as  the  ship's  beam  aft  (see  figure  2).  This  explicitly  assumed  triangular 
ship  geometry  is  the  basis  for  the  computation  of  in  this  research  contribution  and 
the  computer  program. 

A triangular  configuration  provides  a reasonable  approximation  of  the  shape  of  most 
modern  combatants  in  the  critical  locations  near  the  bow  and  stem.  The  small  errors 
that  result  from  using  this  configuration  take  effect  principally  at  small  non-zero  torpedo 
angles  and  to  a lesser  degree  for  near -180  degree  stem  shots.  Typical  beam  Ho -length 
ratios  result  in  an  angle  at  the  bow  of  3-1/2  to  4 degrees  relative  to  the  ship's  centerline. 

A modem  combatant  has  an  angle  of  approximately  10  to  12  degrees.  With  a track 
angle  between  these  values,  the  model  understates  hits  at  the  bow.  Except  at  track 
angles  of  exactly  0,  90,  and  180  degrees,  the  model  overstates  hits  at  the  stem.  These 


-4 


Aim  point: 
center  tine,  amidships 


L = target  ship  length 
B = target  ship  beam 
r = track  angle 
Leff  = effective  length 


Torpedo  track 


FIG.  2:  ASSUMED  SHIP  GEOMETRY 


errors  partially  compensate  for  one  another.  The  model  is  correctly  biased  to  the  extent 
that  torpedoes  traveling  at  a high  relative  velocity  at  angles  of  3 to  12  degrees  relative 
to  ship’s  track  would  be  subject  to  hydrodynamic  forces  at  the  bow  and/or  not  detonate 
because  the  angle  of  incidence  relative  to  hull  plating  is  very  small. 

The  mathematical  approximation  of  for  a triangular  shape,  including  a first- 
order  adjustment  for  relative  motion,  is: 


L x I sin  T I B , u , u \ 

+ I X I V + cos  T I * l 1 + V cos  T) 
5 + COS  T I ' ' 


This  approximation  is  used  in  the  APL  computer  program  described  later. 


Torpedo  Salvos 

To  develop  the  limits  of  integration  for  multiple  hits  on  a target,  the  effective 
length  of  the  target  is  converted  to  an  angular  measure  (the  angle  subtended  from  the 
aim  point  when  the  torpedo  is  at  the  end  of  its  intended  run,  L^/RT)  and  defined 

as  a ratio  to  the  standard  deviation  of  "he  lead  angle  error:  TL  = L^/E^RT  . .The 

coverage  factor,  CF  , is  expressed  as  a multiple  of  TL. 

✓ 

It  will  be  necessary  to  integrate  the  standard  normal  distribution  between  some  limit 
a^  and  b^  , to  compute  the  probability  of  i tracks  intersecting  the  target.  From  this, 

the  probability  of  h hits  can  be  calculated  using  torpedo  reliability. 

The  limits  of  integration,  a.  and  b.  , can  be  represented  algebraic ly  in  terms  of 

TL,  CF,  and  n (number  of  torpedoes  in  the  salvo).  For  convenience  we  will  also  define 
the  spacing  between  torpedoes  as  SP  = (CF  x TL)  / (n-1)  . 

Under  the  assumptions  of  this  model,  a target  is  missed  (no  track  intersections), 
when  the  center  of  the  salvo  is  displaced  from  the  aimpoint  by  at  least  TL/2  + 1/2 
the  salvo  width;  that  is,  (CF  + 1)  x TL/2  = aQ,  » = bQ  (see  figure  3). 


o 8o  bo“°° 


Limits  of  integration 


FIG.  3:  LIMITS  OF  INTEGRATION:  5 TORPEDOES,  NO  HITS 


J3 


We  will  get  one  track  intersection  if  the  center  of  the  salvo  is  displaced  from  target 
center  less  than  a^  = and  more  than  a^  - SP  = (see  figure  4).  Similarly,  we  will 

get  two  track  intersections  if  the  center  of  the  salvo  is  displaced  less  than  h = a , 

but  moire  than  a2  = “ SP  = aQ  - 2 x SP  (see  figure  5).  This  leads  to  the  general  form: 

ai=a()-ix  SP; 

✓ 

bj  = a.^  = aQ  - (i  - 1)  x SP. 

Since  there  are  only  n-1  spaces,  0<isn-l  . However,  these  values  are  not  ex- 
actly correct  because  of  two  limitations: 

First:  If  the  coverage  factor  is  greater  than  1 (i.  e, , the  salvo  spread  is  wider  than 
the  effective  length  of  the  ship),  then  a point  will  come  at  which  a.>0  but  i x SP>TL  . 

That  is,  torpedoes  will  miss  both  ahead  and  behind  the  target  ship  (see  figure  6).  Beyond 
this  point,  the  number  of  track  intersections  alternates  between  i and  i-1  . This 
fluctuation  can  be  visualized  by  thinking  of  the  center  of  the  salvo  moving  from  bow  to 
stem  across  the  target.  At  some  point,  a torpedo  "falls  off'  the  stem  (there  are  now 
i-1  intersections)  and  then  another  "comes  on"  at  the  bow  (i  intersections),  and  so  on. 
After  the  salvo  center  passes  the  target's  mid  point,  the  number  of  intersections 
will  continue  to  fluctuate  for  a while  and  then  diminish  until  eventually  no  tracks  inter- 
sect the  target. 

The  limits  of  integration  are  handled  in  the  following  manner.  When  i x SP>TL  for 
the  first  time,  the  limits  of  integration  for  i track  intersections  are : 

a = a - TL  and  b.  = a.  , , 

i 0 li-l 

After  the  point  when  the  center  of  the  salvo  is  at  a^  - TL  , we  get  one  less  track 

intersection  because  one  torpedo  misses  astern  as  the  salvo  moves  aft.  As  the 
salvo  continues  to  move  aft,  the  next  torpedo  "comes  on"  at  the  bow  when  the  salvo 
center  is  at  aQ  - (i)  x SP  (recall  i x SP>TL) . Thus  integration  between  the  limits 

Sq  - (i)  x SP  and  aQ  - TL  contributes  to  the  probability  of  (i-1)  hits. 

As  the  salvo  moves  aft,  we  get  i hits  until  the  second  torpedo  "falls  off’  the  stem 
when  the  salvo  center  is  at  aQ  - (TL  + SP)  . We  then  get  i-1  hits  until  another 

"comes  on"  at  the  bow  when  salvo  center  is  at  aQ  - (i  + 1)SP  . Continuing  in  this 

pattern,  we  see  that  for  i>(TL-rSP)  , we  end  up  with  alternate  intervals  contributing 
to  the  probability  of  i and  i-1  hits. 


7- 


Two  hits 


a1  *b2 


Limits  of  integration 


Three  hits 


a2=b3 


Limits  of  integration 


Note:  3 x SP>TL.so a3=aQ-TLnotaQ-3SP 


Two  hits 


a3-b4  a2 


Limits  of  int*  ration 


Note:  Mittet  to  both  left  and  right  of  target 


FIG.  6:  INTEGRATION  LIMITS  FOR  CF>1:  6 TORPEDOES,  CF  - 1.8 
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Intervals  of  the  form  aQ  - (i  + j)SP  to  aQ  - (TL  + (j)SP)  contribute  to  the  probabilu| 
of  i-1  hits. 


Intervals  of  the  form  a^  - (TL  + (j)SP)  to  - (i-1  + j)SP  contribute  to  the  proba- 
bility of  i hits. 


That  is 


Prob 

of 

i hits 


s'  - (i-l+0)SP 
N(0, 1) 


■; 


aQ  -(TLHO)SP) 


7 


-(i-l+l)SP 
N(0, 1) 


aQ  -(i-l+2)SP 


aQ  -(TLKl)SP) 


7 


N(0, 1)  +... 


aQ  -(TL+(2)SP) 


and 

Prob  “o:TU<0)SP) 


of  = / 

i-1  J 


aQ  -(TL-Hl)SP) 


hits 


0 


N(0, 1) 
-(i-K))SP 


/ 


N(0, 1) 
-(i+l)SP 


aQ  -(TLrK2)SP 

T N(0, 1)  + ... 

an  -(i+2)SP 


where  N(0, 1)  is  the  standard  unit  normal  distribution.  These  are  still  not  exactly  ca| 
values.  Each  series  must  be  stopped  after  a limited  number  of  terms.  This  leads  to 
second  limitation. 


Second;  We  will  eventually  multiply  the  derived  probabilities  by  2 because  the  aim 
point  was  assumed  to  be  the  center  of  the  target  and  the  standard  normal  distribution 
is  symmetric.  This  means  we  must  integrate  only  over  intervals  to  one  side  or 
the  other  of  target  center  (see  figure  7). 


A check  must  be  made  at  each  step  to  insure  the  limits  of  integration  are  all  non- 
negative. 


The  smallest  integer  value  of  i for  which  i x SPiTL  is  an  upper  limit  on  the 
number  of  intersecting  tracks,  along  with  n (the  number  of  torpedoes  in  the  salvo). 
The  same  constraint  can  also  be  expressed  as  follows:  the  maximum  number  of 
tracks  is  the  smallest  integer  i , such  that  ifen/CF  and  i^n  . 


Reliable  Hits  on  a Target 

Once  the  probabilities  of  0, 1, 2, ...  ,n  intersecting  tracks  are  computed  as  previous! 

described,  the  next  step  is  to  determine  the  probabilities  of  0, 1, 2 n reliable  torpedi 

hits.  To  make  a hit,  a torpedo  must  travel  an  intersecting  track  and  perform  reliably. 

If  we  use  the  symbol  P(T=i)  for  the  probability  that  there  are  i intersecting  tracks  and 
6 for  the  reliability  of  an  individual  torpedo,  then  the  probability  of  h reliable  hits 
P(H=h)  where  h^i  can  be  computed  as: 


P(H> 


,h>=lh 


P(T=i)  x P(H=h  | T=i) 


(5) 
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Assuming  that  die  reliability  of  one  torpedo  is  independent  from  that  of  another,  we 
can  use  a binominal  distribution  for  P(H-h  | T*i) , and  then 


P(H-h) 


If,  «™>|(h)  a-e)<i'h>eh| 


(6) 


Target  Vulnerability 


The  probability  that  a ship  receives  a particular  level  of  damage  Is  the  probability 
of  a particular’  level  of  damage  given  h reliable  hits  times  the  probability  of  h reliable 
hits  summed  over  the  range  of  possible  values  for  h . If  we  symbolize  the  vulnerability 
vector  providing  the  probabilities  of  damage  given  h hits  as  V(h),  h=l,2,...,n  , we 
can  compute  the  probability  of  damage  (PD)  to  a ship  by  the  formula 


PD 


£ 


P(H=h)  x V(h)t 


(7) 


The  damage  probability  can  also  be  interpreted  as  the  expected  fraction  of  attacked  ships 
in  a force  that  are  damaged. 

Simultaneous  Attacks 

To  properly  account  for  simultaneous,  independent  attacks  by  more  than  one  sub- 
marine, the  probabilities  of  reliable  hits  for  each  submarine  must  be  combined  into  a 
probability  of  reliable  hits  by  all  submarines,  before  performing  the  step  indicated  in 
equation  7 above.  For  example,  when  there  are  two  submarines,  we  first  compute 
and  P2(H=h)  for  each  submarine.  Then 


Pj(H=h) 


P(H=h)  = ^ P j(H=i)  x P2(H=h-i) 


(8) 


Although  equation  8 assumes  attacks  by  only  two  submarines,  it  can  be  extended  to  any 
arbitrary  number  of  simultaneous,  independent  attacks. 

COMPUTER  PROGRAM 

The  APL  computer  programs  and  subroutines  in  appendix  A provide  a variety  of 
techniques  for  computing  intersecting  tracks,  hits,  and/or  damage  probabilities  for 
one  or  more  attacking  submarines. 


The  first  program  listed,  and  the  simplest  to  use,  is  Program  ATTACK.  This 
program  is  limited  to  computing  results  for  attacks  by  one  submarine;  it  requires  only 
limited  knowledge  of  the  program  and  APL.  ATTACK  is  a conversational  program  that 
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computes  and  displays  intersecting  track  probabilities,  hit  probabilities,  and/or 
damage  probability.  It  also  displays,  on  request,  the  results  of  such  intermediate 
calculations  as  track  angle,  torpedo  run  time  and  distance,  effective  ship  length,  and 
the  value  of  TL.  Intersecting  tracks  and  hits  are  displayed  as  probability  vectors 
(the  probability  of  0, 1, 2, ...  n tracks  or  hits).  Appendix  B shows  a sample  ATTACK 
run  in  which  user  inputs  are  underlined. 

Appendix^,  also  contains  7 additional  programs  that  do  the  same  computations 
as  program  ATTACK,  but  with  more  flexibility.  These  programs  permit  the  analysis 
of  more  complex  tactical  situations  (i.  e. , multiple  submarines).  These  programs  are 
more  complex  than  ATTACK,  however,  and  familiarity  with  APL  is  essential  to  use 
them  effectively.  The  basic  use  of  each  of  the  programs  is  shown  in  table  1.  Several 
examples  are  shown  in  appendix  C. 
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ATTACK 

✓ 

Computes  intersecting 
tracks,  hits  and 
damage  probability 
for  one  submarine 

Response  to 
conversational 

Inquiry 

Labeled  displays 
aa  requested 

• 

TRACKS 

Computes  probabilities 
of  0,l,2,...,n  inter- 
secting tracks 

A vector  of  target, 
torpedo  and  tactical 
parameters. 

See  page  C-l. 

A vector  of  inter- 
mediate calculations, 
and  intersecting  track 
probabilities . See  page  C-l. 

BITS 

Computes  probabilities 
of  0.1.2. . . ,n  hits 

Track  Intersection 
probabilities  (part 
of  the  output  of 
TRACKS)  plus  torpedo 
reliability 

A vector  of  hit 
probabilities.  See 
page  C-2. 

Elyflgg 

Compute*  probability 

of  damag* 

Output  of  HITS  (or 
SUBS  or  C0M8IHE) 
plus  a vulnerability 
vector  (probability) 
of  damage  given 
0,1,2,...  hits*). 

See  page  C-2. 

The  probability  of 
damage  (expected 
fraction  of  target 
ships  damaged) 

« 

SUBS 

Computes  a track  or 
hit  probability 
vector  for  two  or 
sore  submarines 
firing  from  the 
same  or  from 
symmetric  positions, 
or  for  two  or  more 
salvos  from  the- 
same  (relative) 
position 

The  output  of  HITS 
or  the  probability 
vector  output  of 
TRACKS  and  number 
of  submarines 
(or  salvos). 

The  probability 
vector  for 

0,l,2,...i  hits  or 
tracks  where 
z -number  of  sub- 
marines (or  salvos) 
times  the  size 
of  salvo 

COMBINE 

Combines  twob  prob- 
ability vectors 
from  TRACKS  or  HITS 
into  a single  prob- 
• ability  vector  (to 
combine  hit  prob- 
abilities for  non- 
identical attacks) 

Two  probability 
vectors  of  length 

1^  and  L2 

One  probability 
vector  of  length 

VVl 

* 

BINAP PROX 

Computes  a vector  of 
binomal  probabilities 
(by  normal  approxi- 
mation if  necessary) 

M.P  — tho  binomal 
parameters 

A probability  vector 
of  length  IHl 

NORM 

Computes  areas  under 
the  standard  (unit) 
normal  curve  by 
approximation 

One  or  more  values 
. representing  dis- 
tances from  the  mean 
measured  in  standard 
deviations 

A left-hand 
cuaulative  normal 
probability  for 
each  input 

If  the  vulnerability  vector  docs  not  contain  probabilities  for  all  poaaibla 
numbers  of  hits,  it  ia  extended  by  projecting  a atralght  line  through  the 
laat  two  values  given  until  reaching  an  upper  bound  of  1.0. 

b 

COMBINK  can  be  chained  to  handle  a combination  of  soro  than  two  vectors. 
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V ATTACK  > X » SW . HEAD » R > VUL 
Cl]  ‘WHAT  IS  TARGET  SHIP  LENGTH  IN  FEET?' 

C2]  X«-Q 

C33  'WHAT  IS  TARGET  SHIP  BEAM  IN  fEET?' 

C43  X+-X.0 

C53  'WHAT  IS  TARGET  SHIP  SPEED  IN  KNOTS?' 

C63  X*X»D 

C73  'WHAT  IS  RANGE  TARGET  TO  SUB  IN  YARDS?' 

C83  X«-X.O 

C93  'WHAT  IS  ASPECT  ANGLE  (ANGLE  ON  THE  TARGET  BOW)  IN  DEGREES?' 
C103  X*X.D 

Cl 13  'HOW  MANY  TORPEDOS  IN  A SALVO?' 

C123  X+X'O 

C133  ’WHAT  IS  THE  COVERAGE  FACTOR?  (TORPEDO  SPREAD  + TARGET  SHIP 
LENGTH) ' 

C143  X+X.D 

C153  'WHAT  IS  TORPEDO  SPEED  IN  KNOTS?' 

C163  X*X.O 

Cl-73  'WHAT  IS  TORPEDO  MAXIMUM  RANGE  IN  YARDS?' 

C183  X«-X.D 

C193  'YOUR  ANSWERS  IN  ORDER  WERE: • .ORfXfOR, 'ARE  THEY  OK?' 

C203  -»1XW'N0  ' =5t0 

C213  'DO  YOU  WANT  TRACK  PROBABILITIES  DISPLAYED?' 

C223  SW«-a/'YES  '=5tD 

C233  'DO  YOU  WANT  TRACK  ANGLE.  RUN  DISTANCE.  RUN  TIME.  EFFECTIVE 
LENGTH.  AND  TL  DISPLAYED?'  • 

C243  SW«-SW.a/*yES  '^5t0 

C253  .'DO  YOU  WANT  HITS  DISPLAYED?' 

C263  SW*SW.a/'YES  *=5tD 

C273  'IF  YOU  HAVE  A VULNERABILITY  VECTOR  I"LL  COMPUTE  EXPECTED 
LOSSES. ' 

C283  'INPUT  THE  VECTOR  STARTING  WITH  PROB  OF  LOSS  GIVEN  ONE  HIT.PROB 
OF  LOSS  GIVEN  TWO  HITS.  .-..ETC.  IF  AVAILABLE.  OTHERWISE  HIT  THE 
RETURN. ' 

C293  -»Ox  i 0=+/SW«-SW»  1*pVUL«- » * ' 0 '.lOOtD 
C303  -*Llx\0=+/SWC3.43 

C313  ' I ' ' LL  NEED  TORPEDO  RELIABILITY i ' 

C323  R«-D 

C33D  li:temp«-tracks  X 

C34D  -»N0SH0Txtl=A/(i80. (2.10*7) .0.0. 1 >=6tTEMP 
C353  •*L2x\SWC13=0 

C3A3  'TRACK  PROBABILITIES ' .OR. ( (5  0f0>.10  0» 1 XC&3 ) .OR. 10  5»5*TEMP 
C373  L2:-»L3x\SWC23=0 

C383  'TRACK  ANGLE  RUN  DISTANCE  RUN  TIME  EFFECTIVE  SHIP 

LENGTH'. DR.'  (DEGREES)  (YARDS)  (SECONDS)  (FEET) 

TL(STD.DEV.) ' 

C393  7 1 15  0 15  1 13  0 15  4<iStTEMP 

C403  L31  -»L4x  xSWC33«0 

C413  'HIT  PROBABILITIES' .DR. ((5  0f0>.10  OV I XC63 ) .OR. 10  5tR  HITS  51 
TEMP 

C423  L4!-*0x  iSWC43=,0 

C433  'EXPECTED  LOSS  PROBABILITY!  '*VUL  EL0S8  R HITS  51TEMP 
C443  -»0 

C453  NOSHOT! 'THERE  IS  NO  CHANCE  OF  A HIT  GIVEN  THIS  SET  OF  DATA.' 


V Z«-TRACKS  X$SLrSB>SS>RA» ANrNTf CFr TSrMRfLAfTLtRTIKISPf Hf TIEII IEL 
C13  -»0x\9*pX 

C23  SL«-XC13 

C33  -*N0SH0Tx\0=‘SL4SB<-XC23 

C43  SSOCC33 

C53 4NOSHOTxxO=RA«-XE43 

TT63  AN«-c.XC53*180 

E73  -*NOSHOTx  \0=NT«-XC63 

C83  CF«-XC73 

C93  ->NOSHOTx  \0=TS«-XC83 

C103  -*NOSHOTxto=MR«-XC93 

C113  **NOSHOT  x \ 1 < < SS*TS  > x l oftN 

C123  LA«-“1o(SStTS)  xloftN 

C133  -*NOSHOTx\OiT+-l  + <SS-4-TS>x2oAN+LA 

C143  Z«-180x  (AN4LA)-rOl 

C153  •*L0xtsS=0xE<-*”2.4-0.0006xRA  ’ 

Cl 63  £♦-  ( <E*2)4  ( (0 . 0225X  < lc>AN)  *2)40 « 03045025*  (2oAN)  *2  > + ( <TS+SS>*2)-<  lo 

AN) *2 )*0.5 

C173  L05-»Llx\0=l  I AN-rOl40*RTi-RA*T 
C183  RT«-RAx(lGAN)*lo<c.l)-AN+LA 
C193  L.ll-*NOSHOTx»MR«RT 

C203  EL><T+<SLx|ic.AN4LA>rT«-SBx0.5x|  <SS*TS)42oAN4LA)+14<SS+TS>x2oAN4LA 

C213  Z«-Z»RT * (RT+TSxO. 5626032861 ) »EL»TL«-EL*Ex3xRT 

C223  -*0NEx\<NT*l>vCF=0 

C233  M*(14CF)xTL*2 

C243  SP*<CFxTL)-rNT-l 

C253  -»CFLx  tCFsl 

C263  K<-LTL*SP 

C273  X«-N0RM  M»M-SPx»K 

C283  X*-<1»X)-X»N0RM  M-TL 

C293  I«-0 

C303  AGAIN:XC"1+i>X3*-XC"1+i>X3-(NORM  T*-0rM-SPxK4I«-I4i  >-N0RM  M-TL+ixSP 
C313  -*END*XT=0 

C323  XCpX3*-XCpX34<NORM  T)-NORM  OrM-TL+IxSP 

C333  •♦END*  l0iM-TL4IxSP 

C343  -»AGAIN 

C353  CFLJK«-NT-1+CF=1 

C363  X«-N0RM  M.M-SP*\K 

C373  X*-<lrX)-X»0.5 

C383  ENDtZ*-Z»<NT41>t2*X 

C393  -»0 

C403  0NE*Z«-Z»  < 1-8P)  t < (NT-i)pO) »SP*2x<N0RM  TL*2>-0.5 
C413  -»0 

C423  NOSHOT lZ*<64NT)t 180 t <2p10*7) »0»0»1 

v 
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v Z«-P  HITS  VIN»I 
C13  N«-pV 
C2]  Z4-NtVCl4-i] 

C3D  AGAIN  .Z4-Z+Nt  VC  1 + 1 ]xBINAPPROX  I»P 
C4D  -*AGAINx  iN>  I«-I+l 


v R«-VUL  ELOSS  PROB 
Cl]  •*<Ll»L2»L3)C2+x  < pVUL  >-*PROB] 

C2]  LI  t-*L2x  \ 1==VULCpVUL] 

C3]  -»L3xiO*<>VUL«-<pPROB>  tlLVUL»VULC*>VUL]+<VULCpVUL]-VULC“l+pVUL]>xn. 
PROB 

C4]  L2S  VUL«-(  pPROB)  tVUL  r < pPROB) p 1 
C5]  L3SR<-+/VULxPR0B 


v R4-N  SUBS  P 

Cl]  -*Ox»(N«-LN><1»OpR«-1  0 
C2]  -*Ox  \N=l » OpR«-P 

C3]  R4-P  COMBINE  R 

C4]  •♦3xtl<N<-N-l 


v R«-P  COMBINE  Q 

Cl]  R4-+/C1]  (Or-l  "1  + pP>. MP».  xQ),  <<o,“1>  + pP>  pQ 


v V4-BINAPPROX  YJXJNJPfMJVAR 

Cl]  -»0x\(Pil)vpi0»0pV«-(N+l>t(liP<-YC2])*l»  <NpO>  rl»X*Or  \ N<-LYC  1 ]+0. 5 
C2]  •*3+(^>200)v<N>50)AN>25x(P+l-.P>r  (1-P)+P 

C3]  -»0XV*-(P*X>  x ( X i N>  x ( 1-p) #N-X 

C4]  V«-(  (VARxo2)*“0.5>x*-(  <X-M>*2>+2xVAR<-<l-P)xMi-NxP 

v 


v P 4-NORM  XrN 
Cl]  N4-1,P4-0«*X*,X 

C2]  P*-P»l-(  <4- (2xol) *0.5 >x*~0.5xXCN]*2>x+/0. 31938153  “0.356563782 

1 • 7B1477937  “1.821255978  1 .330274429X (+1+0. 2316419XXCN] ) *15 
C3]  -*2x  i ( pX  > 2N4-N+1 


V 


— 


Note:  Underlines  indicate  user  entries. 


ATTACK 

WHAT  IS  TARGET  SHIP  LENGTH  IN  FEET? 

01 

1000 

WHAT  IS  TARGET  SHIP  BEAM  IN  FEET? 

01 

120 

WHAT  IS  TARGET  SHIP  SPEED  IN  KNOTS? 

01  •" 

20 

WHAT  IS  RANGE  TARGET  TO  SUB  IN  YARDS? 

□ 1 

3500 

WHAT  IS  ASPECT  ANGLE  (ANGLE  ON  THE  TARGET  BOW>  IN  DEGREES? 
□ 1 

135 

HOW  MANY  TORPEDOS  IN  A SALVO? 

□ 1 


6_ 

WHAT  IS  THE  COVERAGE  FACTOR?  (TORPEDO  SPREAD  TARGET  SHIP  LENGTH) 

□ 1 


2*5 

WHAT  IS  TORPEDO  SPEED  IN  KNOTS? 

□ 1 


§0 

WHAT  IS  TORPEDO  MAXIMUM  RANGE  IN  YARDS? 

01 

10000 

YOUR  ANSWERS  IN  ORDER  WERE! 

1000  120  20  3500  135  6 2.5  50  10000 

ARE  THEY  OK? 

YES 

DO  YOU  WANT  TRACK  PROBABILITIES  DISPLAYED? 

YES 

DO  YOU  WANT  TRACK  ANGLE.  RUN  DISTANCE.  RUN  TIME.  EFFECTIVE  LENGTH.  AND 
TL  DISPLAYED? 


YES 

dcT you  want  hits  displayed? 

YES 

IF  YOU  HAVE  A VULNERABILITY  VECTOR  I'LL  COMPUTE  EXPECTED  LOSSES. 

INPUT  THE  VECTOR  STARTING  WITH  PROB  OF  LOSS  GIVEN  ONE  HIT.PROB  OF  LOSS 
GIVEN  TWO  HITS.  ...ETC.  IF  AVAILABLE.  OTHERWISE  HIT  THE  RETURN. 
>5  .75  .1 

I'LL  NEED  TORPEDO  RELIABILITY} 

□ 1 


.85 

TRACK  PROBABILITIES 

0 12  3 

0.20017  0.1 6000  0.63983  0.00000 

TRACK  ANGLE  RUN  DISTANCE  RUN  TIME 
(DEGREES)  (YARDS)  (SECONDS) 

151.4  • 5175  164.0 

HIT  PROBABILITIES 

0 12  3 

0.23857  0.29915  0.46228  0.00000 

EXPECTED  LOSS  PROBABILITY!  0.4962868154 


4 5 6 

0.00000  0.00000  0.00000 
EFFECTIVE  SHIP  LENGTH 
(FEET)  TL(STD.DEV.) 

781  0.7320 


4 5 6 

0.00000  0.00000  0.00000 


APPENDIX  C 

EXAMPLES  OF  APL  PROGRAM  USE 


Example  1 : Torpedo  tracks  (program  TRACKS) 

■ ✓ 

Maximum  torpedo  range  (yards) 

Torpedo  speed  (knots) 

Coverage  factor 

Number  of  torpedoes  in  salvo 

Angle  on  the  target  bow  (degrees) 

Range  to  target  (yards) 

Target  speed  (knots) 

Target  beam  (feet) 

Target  length  (feet) — i 


TRACKS  1000  120  20  3500  135  6 2.5  50  10000 

151.4299402  5175.037571  183.9675558  781.4423676  0.7320379418- 

0.2001704761  0.1599962003  0.6398333226  0000 


Probability  of: 

intersecting  tracks  — I 

intersecting  tracks 

intersecting  tracks 

intersecting  tracks  — 
intersecting  tracks  — 
intersecting  tracks  — 
= n intersecting  tracks 


The  first  values  in  the  output  are: 

Track  angle  (degrees);  run  distance  (yards),  run  time  (seconds), 
. effective  ship  length  (feet),  and  TL,  respectively. 


r 


L 

■ 


Torpedo 

reliabiljjty 


APPENDIX  C 

EXAMPLES  OF  APL  PROGRAM  USE 

Example  2:  Torpedo  hits  and  tracks  (programs  HITS  and  TRACKS) 

— To  remove  intermediate  values  from  TRACKS  output 

Same  as  input  for  TRACKS,  example  1 

I 


.85  HITS  54  TRACKS  1000  120  20  3500  135  6 2.5 

0.2385661559  0.2991542675  0.4622795756  0000 


50  10000 


-1  hit 
-0  hits 


Example  3:  Expected  losses  from  torpedoes  (programs  ELOSS,  HITS,  and  TRACKS) 

Probability  of  damage  given  0 hits 

Probability  of  damage  given  1 hit 

Probability  of  damage  given  2 hits 

I Probability  of  damage  given  3 hits 

Same  as  above,  example  2 


0 .5  .75  1 ELOSS  .85  HITS  54  TRACKS  1000  120  20  3500  135  6 2.5  50  10000 
0.4962868154 


■ Probability  of  damage 


C-2 


APPENDIX  C 

EXAMPLES  OF  APL  PROGRAM  USE 


Example  4:  T orpedo  hits  by  multiple  submarines  (programs  SUBS,  HITS,  and  TRACKS) 

2 SUjpS  .85  HITS  5*  TRACKS  1000  120  20  3500  135  6 2.5  50  10000 


0. 0569 1381074  0.1427361672  0.3100617983  0.2765858156  0.213702406 

00000000 

Probability  of  0,  1,  2 2xn  hits  with  two  submarines  attacking  from  the  same  relative  position. 


Example  5:  Expected  losses  from  a multiple  submarine  attack;  (programs  ELOSS,  SUBS, 
HITS,  and  TRACKS) 

0 0.5  0.75  1 ELOSS  3 SUBS  0.85  HITS  51TRACKS  1000  120  20  3500  135 
6 2.5  50  10000 


0.9251381545 


Probability  of  damage  with  three  submarines  attacking  from  the  same  relative  position. 


Example  6:  Expected  losses  from  a multiple  submarine  attack;  submarines  in  different 

relative  positions  (programs  ELOSS,  SUBS,  HITS,  TRACKS,  and  COMBINE) 

0 0.5  0.75  1 ELOSS  ( 3 SUBS  0.85  HITS  51  TRACKS  1000  120  20  3500  135 
6 2.5  50  10000)  COMBINE  .9  HITS  51  TRACKS  1000  120  20  4000  60  5 1.5  6 

0 8000 


0.9900551757 

Probability  of  damage  with  three  submarines  attacking  from  the  same  relative  position  and  one 
submarine  attacking  from  a different  position  with  a different  torpedo  type. 


